The effect of the spin-orbit interaction on the band gap of half-metals 
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The spin-orbit interaction can cause a nonvanishing density of states (DOS) within the minority- 
spin band gap of half-metals around the Fermi level. We examine the magnitude of the effect in 
Heusler alloys, zinc-blende half metals and diluted magnetic semiconductors, using first-principles 
calculations. We find that the ratio of spin-down to spin-up DOS at the Fermi level can range from 
below 1% (e.g. 0.5% for NiMnSb) over several percents (4.2% for (Ga,Mn)As) to 13% for MnBi. 
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Half-metals are materials showing spin ferromag- 
netism, but with the exotic property of presenting a 
metallic density of states (DOS) only for the one spin 
direction (usually majority), in parallel with a clear 
band gap around the Fermi level Ep for the other di- 
rection. Experimentally, the first of the kind to be re- 
ported were half-metallic Heusler alloys by de Groot et 
al. in 1983>i More recently, other half-metals were found, 
such as Cr02j£ Lao.7Sr .3Mn03j£ or diluted magnetic 
semiconductors. 3 Furthermore, ordered zinc-blende CrAs 
and CrSb were fabricated by molecular beam epitaxy^ 
and relevant calculations suggest these and many similar 
zinc-blende systems to be half-metallici^ So has been 
the case with a variety of Heusler alloys X The rising in- 
terest in half-metallic systems is partly caused by their 
potential applications in spin electronics,- since junctions 
made of half-metals should theoretically present 100% 
magnetoresistance or 100% spin polarized currents. 

Although the ab-initio calculations show the presence 
of half-metallicity in this variety of materials, with a clear 
band gap around Ep for the minority spin, this can be 
partly suppressed by, e.g., defects^ spin excitations at in- 
creased temperature, or non-quasiparticle states^ Even 
at low temperature and in the defect-free case, there is 
still the question of spin-orbit interaction. Although this 
is weak and in many cases can be neglected, in principle 
it should couple the two spin channels so that the pres- 
ence of states around Ep for the majority spin is partly 
reflected in the minority DOS. Thus one cannot have a 
real gap, but a small DOS, its smallness depending on 
the strength of the spin-orbit coupling for each mate- 
rial; this effect could have also implications in the ideas 
about spin-dependent transport and spintronics devices. 
It is the purpose of this paper to give a quantitative pic- 
ture for this effect for a few typical systems, and to ar- 
gue if and how crucially this could affect spin-dependent 
transport. As model systems we have chosen the Heusler 
alloy NiMnSb, the ordered zinc-blende alloys CrAs, CrSb 
and MnBi, and the chemical disorder in diluted magnetic 
semiconductor Gao.95Mno.05As. 

Our first-principles calculations are based on the 
screened Korringa-Kohn-Rostoker (KKR) Green func- 



tion method, incorporating fully relativistic effects (from 
where the spin-orbit interaction emerges). For the de- 
scription of the diluted magnetic semiconductors, the 
coherent potential approximation (CPA) was also em- 
ployed. Details about the method can be found in 
Ref. 0. In the calculations we have used the atomic 
sphere approximation (ASA) , in which it is assumed that 
the potential around each atomic site is spherically sym- 
metric. In the open structures, such as the zinc-blende, 
the large empty interstitial space is accounted for by plac- 
ing at appropriate positions empty spheres (treated as 
atomic spheres with no nuclear charge). In the expansion 
of the Green function in local orbitals a cutoff of Z max = 3 
was taken for the ordered compounds and of l max = 2 
for the CPA calculations. We have employed both the 
fully relativistic (solving the Dirac equation) and the 
scalar relativistic^ methods for comparison. The lat- 
ter includes all other relativistic effects except spin-orbit 
coupling, so that the two spin channels are decoupled 
and half-metals appear with a clean gap for spin-down 
around Ep. 

The DOS n a (E) per spin a is related to the retarded 
Green function G G (f,r';E) via 



n a {E) = -ilm J G a (r,r;E) d 3 



(1) 



Usually for the calculation a very small imaginary part 
it is added to the real energy E and the Green function 
and DOS are evaluated at E + ie. In this way the DOS of 
every eigenstate is no more a delta function of the energy 
but is rather smoothed via a Lorentzian broadening. This 
method has the advantage of needing fewer k points for 
the integration in the Brillouin zone to obtain a smooth 
DOS, saving computational time. Nevertheless there is 
the disadvantage that, in the presence of a band gap, the 
band edges are not sharp and one also obtains a very low 
nonzero DOS within the gap due to the Lorentzian broad- 
ening. In our case this causes a problem, since the DOS 
within the minority gap due to spin-orbit coupling is also 
small and drowned in the background of the Lorentzian 
broadening. To overcome this difficulty one would have 
to go much closer to the real energy axis, which would 
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demand too many k points. We have used an alternative, 
which saves calculation time: Since the Green function is 
analytical in the complex energy plane E z , we are able to 
extrapolate to the real energy axis by knowing the vari- 
ation of the Green function for E z along a line parallel 
to the real axis. To our knowledge, this method was first 
proposed by Gray and Kaplan^ who also report on the 
numerical limitations. The resulting accuracy within the 
gap is much better than that of the traditional method, 
and good enough for our purposes, although sometimes 
numerical fluctuations of the DOS are visible. 

In our calculations we solve the Dirac equation, thus 
treating the spin-orbit term exactly. But for a better un- 
derstanding we give now a brief description of the first or- 
der effect of the interaction in a Schrodinger picture. We 
remind that the spin-orbit coupling of the two spin chan- 
nels is related to the unperturbed potential V(r) around 
each atom via the angular momentum operator L and 
the Pauli spin matrix a: 



1 hldV 



2 r dr 



f SO v so 
\ y so y so 



(2) 



Analyzed in spinor basis, it takes the 2x2 matrix 
form above. We denote the two spin directions with 
| and |, the unperturbed crystal hamiltonian for the 
two spin directions as H°i and H°i, and the unper- 
turbed Bloch eigenfunctions as -. and \& -. Then 

nk nk 

the Schrodinger equation for the perturbed wavefunction 

nk 



reads 

nk nk 



k x t # oi + v& -e)\ 



f =0 (3) 



nk 



The potential terms Vj Q i and VjJ are responsible for flip- 
ping the spin. Within the gap region of the half-metal, 
where no spin-down states of the unperturbed hamilto- 
nian exist, the spin-down solution can be found in first 
order by solving Eq. for (the index (1) stands 

for the first-order solution). The solution reads formally: 



nk 



Hi 



V s i l -E°l 



with E°1 the eigenenergy of the state \E ,U L. We see that 

nk . . . n . k 

in the gap region the spin-down intensity is a weak image 

of the band structure E -, of the spin-up band. Since the 

nk 

DOS is related to l^^il 2 , it is expected that within the 

nk 

gap the DOS has a quadratic dependence on the spin- 
orbit coupling strength: n^E) ~ (VjJ) 2 . But there is 
also a modification by the term — 1/ (Hi + Vj£ — E°V) 
(this is actually the Green function) , which increases the 
weight of the states close to those A; points where the 
unperturbed spin- up and spin-down bands cross, i.e., 



so nk 



,01 



(4) 



E°l — E°j-.. To see this we can rewrite Eq. Q by ex- 

n k n k 

panding the Green function in spectral representation: 
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(5) 



Here, the summation runs only over the band index n' 
and not over the Bloch vectors k', because Bloch func- 
tions with k' ^ k are mutually orthogonal. Close to the 
crossing point E°t — E ^ the denominator becomes small 

nk n'k 

and the bands strongly couple. Then one should also con- 
sider higher orders in the perturbation expansion. Since 
at the gap edges there exist spin-down bands of the un- 
perturbed hamiltonian, this effect can become important 
near the gap edges. 

Now we proceed to the analysis of our results. The 
results for the Heusler alloy NiMnSb, which was found by 
de Groot et al. in 1983 to be half-metallic, 1 are shown in 
Fig-ffl I n t ne top panel, the atom-resolved DOS is shown, 
following a scalar relativistic calculation. The Fermi level 
is within the band gap for minority spin and the material 
is half- metallic. In this scale, the relativistic calculation 
shows no difference in the DOS compared to the scalar 
relativistic. In the middle panel the gap region is shown 
in a relativistic calculation. Again, no minority DOS is 
seen in this scale. But when we magnify the minority 
DOS in the lower panel, it becomes clear that the DOS 
is nonzero in the gap region — it is of the order of 0.5% of 
the majority DOS in the vicinity of Ep. We have checked 
the DOS in this scale for the scalar relativistic calculation 
and have found it to be zero, thus we conclude that we see 
here a spin-orbit effect and not an artificial Lorentzian 
broadening due to the imaginary part of the energy in 
the calculation. The fluctuations of the Mn DOS (full 
line) are due to numerical noise of the k integration and 
extrapolation scheme, and mark its limit of accuracy. 

The next example is the ordered half-metallic alloy 
CrAs in the zinc-blende structure. This was invented 
by Akinaga and co-workers 4 in 2000 and found to be 
half- metallic. Our results arc presented in Fig. [3 In the 
top panel, the scalar relativistic atom-resolved DOS is 
shown. The Fermi level is within the minority-spin gap, 
close to the conduction band edge. (The latter is mainly 
of d character, containing the wavefunctions of the e g 
subspace 6 ). Focusing in the range of the gap, we can see 
the DOS in more detail in the middle panel in a relativis- 
tic calculation. Below —1 eV one recognises the majority 
d states of e g character^ Around —0.9 eV the DOS has 
a minimum, as the e g states end and the d states of tig 
character start, forming wide bands. In this scale no 
spin-orbit effect is visible in the minority gap. But if we 
focus further in energy (lower panel) it is clear that the 
structure of the majority DOS is reflected in the minority 
DOS: one can recognize the peak of the e g states below 
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FIG. 1: DOS of NiMnSb. Top: scalar relativistic (I. refers 
to interstitial volume of the empty sphere). Middle: fully 
relativistic, focused in the minority-gap region around Ef- 
In this scale, the spin-orbit effect is not visible. Bottom: fully 
relativistic, focused on the minority states within the "gap". 
The fluctuations reflect the numerical accuracy of the method 
(see text). Around Ef, the spin-down DOS is of the order of 
0.5% of the spin-up DOS. Close to the "gap" edges, the spin- 
orbit induced minority DOS increases strongly as explained 
in the text. Negative numbers in the DOS axis correspind to 
minority spin. 



— 1 eV, then the minimum and finally the presence of the 
t2 g band above —0.8 eV. This is the effect of the term 
VjJ g in Eq. J3J . But we also observe the rapid increase 
of the DOS as the energy approaches the band edges; this 
comes from the term -l/(H l + V s l Q l - E°l). The minor- 
ity DOS in the middle of the gap is again small, giving 
a minority/majority DOS ratio of 0.2%. In a similar cal- 
culation for zinc-blende CrSb (first fabricated by Zhao 
et alm&) we find that the minority DOS in the middle of 
the gap is around 0.7% of its majority counterpart. This 
rise is expected, as the spin-orbit interaction is stronger 
for the hp states of Sb than for the Ap of As. (The trend 
can be seen for example in a systematic study of the 



FIG. 2: DOS of zinc-blende CrAs. Top: scalar relativistic. 
Middle: fully relativistic, focused in the minority-gap region 
around Ef- In this scale, the spin-orbit effect is not visi- 
ble. Bottom: fully relativistic, focused on the minority states 
within the "gap". The reflection of spin-up into spin-down 
DOS is clear at the peaks at —1 eV and —0.8 eV. In the mid- 
dle of the gap, the spin-down DOS is of the order of 0.2% 
of the spin-up DOS. Negative numbers in the DOS axis cor- 
respind to minority spin. 



spin-orbit cross section of sp impurities in Mg in Ref. 0, 
where Sb is found to have four times higher spin-orbit 
cross section than As. Since the spin-orbit interaction is 
in first approximation an atomic property, the trend is 
similar here.) 

A recent theoretical investigation^ has suggested the 
ordered zinc-blende alloy MnBi as a good candidate for 
half-metallicity. Although this has not yet been fabri- 
cated, we take it as a good example to show the spin-orbit 
effect, since Bi is a heavy element with strong spin-orbit 
interaction because of the 6p states. The scalar and fully 
relativistic DOS are shown in Fig. [3] It is evident that 
here the spin-orbit interaction has a much stronger influ- 
ence. The minority DOS at Ep is of the order of 13% of 
its majority counterpart, a huge quantity in comparison 
to the analogous percentage in CrAs or NiMnSb. The 
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FIG. 3: DOS of zinc-blende MnBi. Top: scalar relativistic. 
Middle: Fully relativistic. Bottom: fully relativistic, focused 
on the states within the minority "gap". At Ef, the spin- 
down DOS is of the order of 13% of the spin-up DOS. This 
higher value is expected since the 6p valence states of Bi result 
in a strong spin-orbit effect. Negative numbers in the DOS 
axis correspind to minority spin. 



relativistic effect shows up also in the orbital magnetic 
moment, being 0.11/is and mostly concentrated at the 
Mn atom, while the spin moment remains at 4.0/Lts per 
unit cell as in the scalar and non-relativistic case. We 
note that, although the spin-orbit interaction is strong 
here, the majority spin is still dominating at Ep. Thus 
MnBi would be a useful material for spin electronics de- 
spite the spin-orbit effect, as its spin polarization at Ep 
is robust against moderate volume change as pointed out 
in Ref.H 

As a final example we consider the diluted magnetic 
semiconductor (DMS) Gao.95Mno.05As. This material, 
together with other similar DMS, has been extensively 
studied in recent years, because it presents ferromag- 
netism at such low Mn concentrations (discovered by 
Ohnc 3 in 1998). Ideally, the Mn impurities substitute 



Ga atoms in the GaAs lattice. Then the Mn 3d states 
of the t2 g representation (d xy , d yz , and d xz ) hybridize 
with the p states of the As neighbors, forming bond- 
ing and antibonding p-d hybrids and from them states 
within the GaAs gap for the spin-up direction. Thus 
the alloy becomes half-metallic with Ep close to the va- 
lence band edge. Although there is no consensus yet on 
the physical origin of the ferromagnetic ordering, which 
seems to depend on the energetic position of the Mn d 
statesjii we have proceeded using the LDA result. We 
find that the spin-down DOS at Ep is 4.2% of its spin- 
up counterpart if we include spin-orbit effects, while the 
system is half-metallic in the nonrelativistic solution. To 
understand this relatively high percentage of spin-down 
states, we observe that Ep is very close to the valence 
band edge. There, the p states of As give an impor- 
tant spin-orbit effect. This is known to split the valence 
band edge of GaAs in two subspaces, one with total (spin 
+ orbital) angular momentum j = 3/2 and one with 
j = 1/2, so that the valence band edge has no pure spin 
character i& Similar calculations on Mn-doped ferromag- 
netic zinc-blende GaN gave us no spin-orbit effect at Ep 
within numerical accuracy. The difference to (Ga,Mn)As 
is that the 2p valence states of N compared to the Ap of 
As (i) have a much lower spin-orbit interaction, and (ii) 
are lower in energy so that the Mn impurity d band is 
fully within the band gap, away from the N-dominatcd 
valence band states, and the spin-orbit effect weakens. 

In summary, we have performed first-principles calcu- 
lations in order to investigate the effect of the spin-orbit 
coupling to the spin-down band gap of half- metallic sys- 
tems. As typical half-metallic systems, describable by 
the LDA, we have chosen the Heusler alloy NiMnSb, 
the ordered zinc-blende alloys CrAs, CrSb, and MnBi, 
and the diluted magnetic semiconductors (Ga,Mn)As and 
(Ga,Mn)N. We find that the majority-spin states are 
partly reflected into the minority band gap. The intensity 
of the DOS for minority electrons in this energy region 
depends mainly on the strength of the spin-orbit cou- 
pling, heavier sp elements resulting in a higher ratio of 
minority/majority spin DOS. Thus we see a trend of in- 
creasing min./maj. DOS ratio in the middle of the gap 
as we go from CrAs (0.2%) to CrSb (0.7%) and finally 
to MnBi (13%). In NiMnSb we find a ratio of 0.2%. 
Also, it is important how deep the Fermi level lies within 
the gap, since majority states close to the gap edges are 
more drastically spin-flipped than states deep in the gap. 
This results in high minority-spin DOS (4.2% ratio) for 
(Ga,Mn)As where Ep is very close to the valence band 
edge. We have also explained these effects within a first- 
order approximation. For transport applications in spin- 
tronics we conclude that even in compounds with strong 
spin-orbit effect (MnBi or (Ga,Mn)As) the minority DOS 
at Ep is still dominated by its majority counterpart (even 
in MnBi the ratio is only 13%), whence current should 
be transported practically only by majority states. Thus 
in applications other effects, such as gap states due to 
impurities^ stalking faults, or interfaces, could be more 
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important to eliminate by sample improvements. 
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